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ABSTRACT 



Aims. We focus on an analysis of galaxies of high stellar mass and low metallicity. Evidence of recent merger events in their optical 
images allow us to classify galaxies into either disturbed or undisturbed, and to study galaxy properties such as morphology, colours, 
stellar populations, and global environment for the different metallicity ranges and disturbance classes. 

Methods. We cross-correlated the Millenium Galaxy Catalogue (MGC) and the Sloan Digital Sky Survey (SDSS) galaxy catalogue 
to provide a sample of MGC objects with high resolution imaging and both spectroscopic and photometric information available in 
the SDSS database. For each galaxy in our sample, we conducted a systematic morphological analysis by visual inspection of MGC 
images using their luminosity contours. The galaxies are classified as either disturbed or undisturbed objects. We divide the sample 
into three metallicity regions, within wich we compare the properties of disturbed and undisturbed objects. 

Results. We find that the fraction of galaxies that are strongly disturbed, indicative of being merger remnants, is higher when lower 
metallicity objects are considered. The three bins analysed consist of approximatively 15%, 20%, and 50% disturbed galaxies (for 
high, medium, and low metallicity, respectively). Moreover, the ratio of the disturbed to undisturbed relative distributions of the 
population age indicator, D„(4000), in the low metallicity bin, indicates that the disturbed objects have substantially younger stellar 
populations than their undisturbed counterparts. In addition, we find that an analysis of colour distributions provides similar results, 
showing that low metallicity galaxies with a disturbed morphology are bluer than those that are undisturbed. The bluer colours and 
younger populations of the low metallicity, morphologically disturbed objects suggest that they have experienced a recent merger with 
an associated enhanced star formation rate. 

Key words, galaxies: formation - galaxies: evolution - galaxies: abundances - galaxies: interactions. 



1. Introduction 

In the most commonly accepted cosmological paradigm, galaxy 
interactions and mergers play a crucial role in determining 
galaxy properties and are considered as one of the main mecha- 
nisms by which galaxies experience significant changes in mor- 
phology, stellar population content, and star formation activ- 
ity (e .g., iBarton et al.1 120001: iLambas et al.l I200I lAlonso eTaP 
I2OO6I) . Thus, it is generally accepted that the different merger 
histories of galaxies define their evolution, origin, and present 
day properties. 

I 

It is also important to consider the chemical properties 
of galaxies. These can provide fossil records of the ir his- 
tory of formation (iFreeman & Bland-Hawthonil l2002h since 
the metallicity content of a galaxy is expected to depend 
strongly on its evolutionary state, namely, how and when was 
the gas transformed into stars. The relation between interac- 
tions/mergers and chemical properties have been studied by dif- 
ferent au thors ( Donzelli & Pastoriza"2000*: 'Marauez et a l.ll2002t 
Fabbian o~t al.ll200 4; Kewley et al. 2006a; Michel-Dans ac et al.l 
,2008). The underlying idea in these studies has been that a 
close companion can induce gas inflows w hich lower the metal- 
licity in the central region s of galaxies (iKewlev et al.l l200"6al: 
iMichel-Dansac et al.ll2008l) . 



A stellar mass-gas metallicity relation ha s been well estab- 
lished (hereafter MZR, iTremonti et al.ll2004 . Although it is a 
clearly observed trend of increasing gas-phase metallicity with 
stellar mass, there is considerable scatter in the relation, which 
could be attributed to the particular star formation history of 
galaxies. Accepted theories for the origin of the MZR have as 
a central proposition that efficient galactic outflows remove met- 
als from galaxie s with shallow potentia l wells (Larson 
lDalcanton ll2007tlFinlator & Davell2008l) . In this context, merg- 
ers and interactions could play an important role in determining 
the shape and scatter of the MZR. 



iKewlev et al.l (l2006ah detected a shift in the luminosity - 
metallicity (LZ) relation towards lower metallicities by » 0.2 
dex for galaxy pairs of a given luminosity, compared to a control 
sample. The spectra analysed in this work corresponded only to 
the central 10% of the galaxy and the authors interpreted this 
result as a signature of metal-poor gas being funneled into the 
center of the galaxies. Following this approach, Ellison et al. 
(2008) studied the LZ relation of 1716 galaxies with compan- 
ions with radial velocity AV < 500 km.s"' and projected sepa- 
ration rp < 80 kpc h and they found an offset to lower metal- 
licities (by » 0.1 dex) for a given luminosity in pair galaxies. 
iMichel-Dansac etal] (l2008b studied the stellar mass - gas metal- 
licity relation of galaxies in close pairs (morphologically classi- 
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fied according to the strength of the interaction signatures) and in 
isolation taken from the Sloan Digital Sky Survey Data Release 
4 (SDSS-DR4). The authors measured differences in the metal- 
licity of galaxies at a given stellar mass, between pairs showing 
signs of strong interactions and galaxies in isolation. Thus, the 
mass-metallicity relation differs between systems with a com- 
panion and galaxies in isolation. Interacting galaxy pairs of high 
stellar mass are found to have systematically lower metallicity 
values than the mean MZ relation, which implies tha t gas inflows 
have b een induced by the interaction. In addition, IPeeples et al.l 
(|2009|) studied 42 outliers (massive, metal-poor galaxies) of the 
MZR in the SDSS, finding that these objects are in general ex- 
tremely blue and star-forming, and show signs of morphological 
disturbances. 

On larger scales, the MZR appears to be only a wea k func- 
tion of galaxy environment (Mouhcine et al. 2007; Coop er et al.l 
12008 ). The gas-phase oxygen abundance of galaxies (at a given 
stellar mass) in the richest environments is only ~ 0.05 dex 
higher than those in the poorest environments, and these trends 
might account for about 15% of the scatter in the MZR. 
[Ellison et al. (2009) explored the impact of cluster membership 
and local-density on the MZR, confirming that galaxies in clus- 
ters are slightly more metal-rich that in the field but that this 
effect is driven by local overdensity and not simply cluster mem- 
bership. 

Given that the majority of galaxies lie close to the best-fit 
mass-metallicity relation, we may be able to learn about the gas- 
phase metallicity evolution of galaxies by studying the prop- 
erties of galaxies that do not follow this relation. Therefore, 
joint statistical analysis of mergers/interactions and chemical 
properties could help us to shed light on the efficiency and ex- 
tent of the processes mentioned above. To this aim, we con- 
sider star-forming, high-stellar-mass galaxies that have in gen- 
eral a well-defined morphology in contrast to smaller systems. 
It is therefore much easier to identify a recent merger event 
in these objects. We study high resolution B-band images of 
high mass (10 < log M./Mq < 10.8) galaxies ext racted from 
the Millennium Galaxy Catalogue (hereafter MGC, iLiske et al] 
[2003; Driver et al. 2005). We classify the galaxies as either dis- 
turbed or undisturbed and we explore statistically the effects of 
recent merger events on oxygen abundance, colour, and stellar 
age indicators derived from SDSS. 



2. Data 

Our data consist of galaxi es from the Millennium Galaxy 
Catalogue dLiske et al.ll2003l; [Driver et al."2005') and th e Sloan 
Digit al Sky Server Data Release 4 ( Adelman-McCarth v et alj 
I2006h . 

The Millennium Galaxy Catalogue is a wide medium-deep 
B-band imaging survey along the celestial equator that was ob- 
tained using the Wide Field Camera on the 2.5-m Isaac Newton 
Telescope on La Palma. The survey covers 37.5 square degrees 
with Z?-band magnitudes in the range 13 < B < 24. The sur- 
vey region coincides with both the 2dF Galaxy Redshift Survey 
northern strip (2dFGRS, Colless et al. 2001) and the Sloan 
Digital Sky Survey region ( Abaz aiian et al.ll2003l) . 

The SDSS is a photometric and spectroscopic survey that 
will cover approximately one-quarter of the celestial sphere and 
collect spectra for more than one million objects. The imaging 
portion of the fourth release SDSS comprises 6670 square de- 
grees of sky imaged in five wavebands (m, g, r, i, and z) contain- 
ing photometric parameters of 180 million objects. The spectra 




log{M.) 

Fig. 1. Stellar mass - gas metallicity relation of a selection of the 
SDSS galaxies (see text). The vertical dotted lines delimit the 
mass range that we study in detail in this paper. The blue points 
correspond to our sample of MGC-SDSS galaxies. 



are obtained from 3 arcsec diameter fibers projected on the sky 
and the coverage is 3800-9200 A. 

The SDSS value-added catalog provides several spectro- 
scopic parameters of physical relevance. However, SDSS imag- 
ing is relatively shallow and often taken in poor seeing condi- 
tions. Therefore, because we wish to classify galaxies accord- 
ing to their degree of disturbed morphology we use a subset of 
SDSS galaxies in common with MGC, which provides deeper 
and higher resolution imaging data. The quality of the images 
in MGC is a major asset for a detailed study of galaxy mor- 
pho logy. Th e median seeing of the MGC survey, 1.3 arcseconds 
(Lis ke etaP 12003 ). combined with its deep surface brightness 
limit of fiB = 26 mag arcsec"^, provides much higher quality 
and more clearly resolved images than SDSS images of surface 
brightness limit 20 mag arcsec"^, which were obtained with a 
median seeing of 1.5 arcseconds. 

We cross-con-elated the MGC and SDSS DR4 galaxy cat- 
alogs to obtain a sampl e of MGC objects w ith oxygen abun- 
dances 12 + log(0 /H) (Tremonti et al."2004'), stellar mass M„ 
and 4000-A break strength D„(4000) (Kauffmann et al. 2003t), as 
well as basic information in SDSS data (e.g., redshifts, concen- 
tration parameter). 

We characterized the local environment of galaxies by defin- 
ing a projected density parameter, Z5. This parameter is calcu- 
lated by using the projected distance to the 5''' nearest neighbour. 
Eg = 5/(ndf). Neighbours were chosen to have luminosities 
brighter than Mr < -19.5 and radial velocity differences lower 
than 1000 km.s"' (Balogh et al. 2004). Kauffmann et al. (2004 
estimated the local density by counting galaxies within cylin- 
ders of 2 Mpc in projected radius and + 50 km.s ' in depth, in 
a complete sample of galaxies from SDSS. [Alonso et aL I (l2006l) 
found a good correlation signal between both density estimators 
indicating that either of these two definitions are adequate for 
characterizing the environment of galaxies. Studies of the rela- 
tionship between gas-phase oxygen abundance and environment 
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have used simi lar parameters, for i nstance the S3 projected den- 
sity parameter (ICooper et alJl2008 '). or a density estimator based 
on the average of the projected distan ces to the fourth and fi fth 
nearest neighbour within 1000 km. s"^ dMouhcine et al.li2007h . 

We imposed an upper redshift cutoff at z < 0.1 to ensure 
both high completeness for both the SDSS and MGC data, and 
sufficiently high angular resolution. Additionally, we considered 
a minimum redshift z > 0.04 to avoid a too low spatial coverage 
by the SDSS fibers. Thus, for the range adopted (0.04 < z < 0.1) 
the 3 arcsec SDSS fiber probes approximatively the central 2.5 
kpc to 6 kpc of the galaxies. 

We restricted the present analysis to high stellar mass ob- 
jects in the range 10 < log M./M© < 10.8. The high stellar 
mass range adopted is consistent with rather extended galaxies 
so that metallicity measures do not include the external parts 
of the galaxies. Moreover, the rather narrow stellar mass range 
adopted ensures that our subsample of galaxies is homogeneous 
and adequate for the detection of substructures. 

To check for any potential AGN contamination in our sam- 
ple, we analysed the BPT diagnostic diagrams [Oiii]/Hy8 ver- 
sus (vs.) [Nii]/Hq', and both [Oiii]/H ^ vs. [OiJ /Hg and 
[Oiii]/H^ vs. [S II] /Ha proposed bv lBaldwin et al.l (il98ll), in 
all case s using the classification criteria given in Kewlev et al.l 
(I2006bl) . We found that a small percentage of the objects have 
line ratios indicative of their spectra being possibly associated 
with composite AGN/H 11 regions. To avoid this possibility, we 
removed these ambiguous cases from the analysis. 

A possible concern with SDSS spectroscopy is the signal-to- 
noise ratio of the relevant lines associated with gas-phase oxy- 
gen abundance determinations. We checked that in our sample 
all galaxies have a signal-to-noise ratio (S/N) of at least 5 in the 
lines Ha, H/3, and [N Ii]/16584, and at le ast 3 in the [O Ii i1^5007 



line following the selection criteria of 'Trem onti et alJ (l2004t) . 
Thus, despite metallicity determination being a complex task 
with a rather large associated uncertainty, our analysis is not 
likely to be affected by systematic biases (see also discussion 
in Sect.[X2]l. 

With the above restrictions, our final sample comprises 191 
high stellar mass galaxies. In Fig. [1] we plot the MZR for the 
total sample of SDSS galaxies in the redshift range 0.04 < z < 
0.1 satisfying the selection criteria adopted (e.g., AGN removal, 
S/N), and our sample of high stellar mass MGC-SDSS galaxies. 
Figure|2]shows the 12 + log(0 /H) distribution in both catalogs 
in 3 mass bins within the mass range selected previously. It can 
be appreciated that the MZR of our MGC-SDSS sample is truly 
representative of a MZR determined for the whole SDSS in the 
selected mass range. 



3. Analysis 

3.1. Visual classification of galaxies into disturbed and 
undisturbed objects 

For each galaxy in our sample we conducted a systematic mor- 
phological analysis by visual inspecting MGC images and study- 
ing the luminosity contours obtained with IRAF routines applied 
to the reduced images provided in the catalogue. Galaxies are 
classified as: 

a) Disturbed: these objects show externally-triggered distor- 
tions indicative of a recent strong tidal interaction or merger. 
Some of these objects exhibit arcs, shells, ripples, tidal de- 
bris, warps, rings, extremely asymmetric light distributions, tidal 
tails, multiple nuclei, and strong substructure inside a common 
body. We note that these galaxies are not two interacting systems 



(1 0.0 <log(M/MO)< 10.27) 
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Fig. 2. Oxygen abundance distributions: dashed lines correspond 
to the SDSS Tremonti et al. data, and solid lines to our sample 
of MGC-SDSS galaxies in three stellar mass bins. 



in a ongoing merger, but a single galaxy with signs of being a 
merger remnant. 

b) Undisturbed: objects with smooth luminosity contours. 

To assess the advantage of using MGC images for the present 
analysis, we compare in Fig.[3]Z?-band MGC and g-band SDSS 
images for a sample of disturbed galaxies where it is clearly ap- 
parent that there is higher reliability when assigning a disturbed 
morphology in the higher quality MGC data. It is clear that sev- 
eral cases would not be classified as disturbed in SDSS, so that 
the improvement provided by using MGC data is evident. We 
also used SDSS data to analyse the close local environment of 
each galaxy, since MGC images have a too small field of view. 
Since we are interested in analysing the effects of past merger 
events, we excluded a few potentially interacting pairs in the 
sample. 

We considered the uncertainty in the visual classification of 
galaxy morphology by independent determination by different 
observers. By cross-checking the results, we estimate a 95% co- 
incidence indicating the reliability of our classification scheme. 

Visual classification is by nature a subjective task. To as- 
sess the reliability o f our cl assification, we derived the asym- 
metry A parameter dConselice et al. 2000) for the MGC im- 
ages Different authors (e.g.. lDe Propris et al.ll2007l: I Jo gee et al] 
120091) have shown that pairs and galaxy mergers exhibit larger 
asymmetries. In our sa mple, for each galaxy, we use SExtractor 
(iBertin & Arnoutsl 'l996) to calculate the A parameter, finding 
that galaxies classified as disturbed by visual inspection have 
a larger asymmetry parameter (A ^ 0.30), and undisturbed ob- 
jects have lower asymmetry values (A ^ 0.30). More precisely, 
we find that 72% of disturbed galaxies exhibit A <: 0.3, whereas 
all the galaxies classified as undisturbed have A ^0.3. 

For the stellar mass range explored, we found 46 strongly 
morphologically disturbed galaxies, 123 undisturbed galaxies, 
and 22 galaxies with a close companion. The redshift distribu- 
tions of disturbed and undisturbed galaxies is shown in Fig. |4] 
The great similarity between these distributions is an indication 
that our classification is not affected by angular resolution. 




Fig. 3. Images and surface brightness contours of typical massive disturbed galaxies in both MGC (left columns) and SDSS (right 
columns) catalogs. The image sizes are 33 x 33 arcsec^. 



3.2. Oxygen abundances, colours, and star formation 
histories of disturbed and undisturbed galaxies 

We analysed the oxygen abundance and stellar mass in our sam- 
ple of high mass MGC-SDSS galaxies. We show the MZR of 
these galaxies in Fig.|5] We considered two threshold curves dis- 
placed with respect to the Tremonti et al. fit, which separate the 
sample into three regions with low, medium, and high metallicity 
defined by: 

high: 12 + log(0/H) > -1.492 + 1.847(logM,) - 
0.08026(logM,)2-H0.078; 

low: 12 + log(0/H) < -1.492 + 1.847(logM.) - 
0.08026(logM.)^-0.03; 

and a medium metallicity range between these two. Adopted 
thresholds divide the sample into approximatively 1/4, 2/4, and 
1/4 the total sample of galaxies (representing the high, medium, 
and low metallicity ranges respectively). 

As shown in Fig. |5] and Table[Tl the relative fractions of dis- 
turbed galaxies are statistically significantly different in the three 



metallicity bins explored. This remarkable tendency for the rel- 
ative number of disturbed galaxies to increase with decreasing 
metallicity is also evident from visual inspection of Fig. |5] To 
quantify these results, we also calculated the percentages of dis- 
turbed galaxies for the three different metallicity bins defined 
previously. At medium and high metallicities, we find that the 
percentages of disturbed objects are not very diff'erent (14% and 
21% respectively, cf Table[T]i. However, at low metallicities, the 
disturbed galaxy fraction increases remarkably (52%). 

To check for possible errors in the metallicities of galax- 
ies in our samp le, we recalculated 12-1- log(0/H) using the 
iPettini & Pageir(l2004 method, where the oxygen abundance is 
given by 12 + log(0 /H) = 9.37 + 2.03 x N2 -h 1.26 x Nl^ + 
0.32 X N2^ and N2 = log([N Ii]i6584/H a). This method is par- 
ticularly suited to our sample since all galaxies have a S/N at 
least 8 in the emission lines used to compute N2. Moreover, 
the N2 range of the galaxies studied (-0.7 < N2 < -0.29) is 
within the domain of validity for this equation. We considered 
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Table 1. Number of galaxies classified as disturbed galaxies, undisturbed galaxies, and galaxies with a close companion pairs, in our 
three ranges of metallicity, and percentages (and respectives sampling errors) of disturbed galaxies with respect to the total number 
of disturbed and undisturbed galaxies 





disturbed 


undisturbed 


with close comp. 


total 


% disturbed 


total 


46 


123 


22 


191 


27.2 ± 3.4 


high 0/H 


6 


35 


6 


47 


14.6 ±5.8 


medium O/H 


18 


68 


13 


99 


20.9 + 4.4 


lowO/H 


22 


20 


3 


45 


52.3 ±7.7 




0.02 0.04 0.06 0.08 0.1 0.12 10 10.5 1 1 

z log(M./M„) 



Fig. 4. Redshift distribution of disturbed (solid line) and undis- 
turbed galaxies (dashed line). 



the three metallicity regions defined previously but shifted by 
-0.35 dex since this m ethod tends to infer low er values of metal- 
licities than those of Tremonti et al] (l2004l) (see for instance 
iKewlev & Eliisoi3l2008ir With these new estimates, we obtain 
similar numbers of galaxies in each metallicity bins, as shown 
in Fig.|6] We also used a subsample of 82 galaxies with data of 
S/N of at least 8 in the strong emiss ion lines and oxygen abun- 
dances measured bv lTremonti et al.l ( 12004.) . finding very similar 
percentages (63, 22, and 17%) of disturbed galaxies in the low, 
medium, and high metallicity bins. Then, this trend of increasing 
the disturbed galaxy fraction at lower metallicity appears robust 
and is not likely to be caused by errors in the determination of 
oxygen abundances for these galaxies. 

Studies of galaxy properties and their dependence on envi- 
ronment are important for understanding the role of mechanisms 
driving the evolution of galaxies. We therefore explored the pos- 
sible influence of the environment of each galaxy in our sample 
on the results. To this aim, we analysed the local density param- 
eter S5, defined previously, for all galaxies in our sample. The 
results of this analysis are shown in Fig. |7] where it can be ap- 
preciated that disturbed and undisturbed galaxies populate simi- 
lar environments indicating that the effects are not caused by the 
particular location of these systems. Thus, we conclude that the 
strong trend of increasing percentages of disturbed galaxies as 
a function of decreasing metallicity is not likely to be biased by 
environmental effects. 



Fig. 5. Distribution of 12 -H log(0/H) and log(M,/M0) values 
for galaxies with stellar masses in the range log(M,/MQ) ^ 10.0 
to 10.8. The circles correspond to disturbed galaxies and trian- 
gles to undisturbed galaxies. The solid and dashed lines indicate 
the limits of the high, medium, and low metallicity ranges. The 
solid line corresponds to the Tremonti et al. (2004) fit to the MZ 
relation shifted by -hO.078. The dashed line corresponds to this 
fit shifted by -0.03 dex. 



For completeness in our analysis, we also studied the light 
concentration parameter (i.e., the ratio of Petrosian 90% to 50% 
r-band light radii, C = ?"9o/'"5o) of the galaxies in the sample 
to search for possible differences in our sample between dis- 
turbed and undisturbed objects. The concentration index value 
of C = 2.5 was adopted to segregate concentrated, bulge-like 
(C > 2.5) galaxies from more extended, disc-like (C < 2.5) sys- 
tems. In Fig. |7] the similarity between the C parameters across 
the samples is clearly evident, showing that the observed dis- 
turbed to undisturbed ratio dependence on metallicity is prob- 
ably not associated with systematic differences in the bulge-to- 
disc ratios of the galaxies. To quantify the significance of the 
differences between the distributions in Fig.|7] we calculated the 
occurrence of C < 2.5 in galaxies with low metallicities in 20 
randomized samples. In each sample, the C parameter was re- 
assigned randomly, in disturbed and undisturbed objects of low 
metallicity. Since we found that in 6 out of 20 randomized sam- 
ples the signal is as strong as in the data, we conclude that no 
statistically significant tendency associated with disturbed and 
undisturbed galaxies can be identified for the concentration pa- 
rameter. 
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C = (r90/r50) 



Fig. 6. Distribution of 12 + log{0/H) calculated through the 
Pettini & Pagel (2004) method, and log(M./Mo) values for 
galaxies with stellar masses in the range log(M,/MQ) x 10.0 
to 10.8. Symbols are the same as in Fig.|5] The solid and dashed 
lines indicate the limits of the high, medium, and low metallicity 
ranges. With these abundances, the solid line corresponds to the 
Tremonti et al. (2004) fit to the MZ relation shifted by -0.272. 
The dashed Une corresponds to this fit shifted by -0.38 dex. 



That the relative fraction of disturbed to undisturbed galax- 
ies increases for low metallicity objects provides important evi- 
dence that a merger process induces low metallicity gas inflows. 
This phenomenon should also be reflected in the mean stellar 
age population of disturbed objects as a f unction of meta lUc- 
ity. The D„(4000) spectral index defined bv lBruzuall (Il983l) is a 
suitable indicator of the current luminosity-weighted mean stel- 
lar age since young stellar populations produce almost no metal 
absorption just shortward of 4000A. Therefore, galaxies dom- 
inated by recent episodes of star formation exhibit a small in- 
dex, while the spectra of galaxies dominated by older popu- 
lations exhibit strong metal lines in absorption and thus larger 
D„(4000) values. We analysed the D„(4000) population age in- 
dicator in the 3 gas-phase oxygen abundance ranges. The results 
are shown in Fig. [8] where we plot the fraction of D„(4000) 
for disturbed objects normalized to that for undisturbed galax- 
ies (fraction{D„(4000)(disturbed)fD„(4000)(undisturbed))).We 
find that at both medium and high metallicity, both disturbed and 
undisturbed objects have a similar D„(4000), with a slight ten- 
dency for disturbed objects to have older stellar populations. It 
can also be appreciated that low O/H values objects, where the 
disturbed population is important, have substantially younger 
stellar populations than their undisturbed counterparts. Thus, 
our morphological analysis indicates that in high stellar mass 
galaxies, low metallicity values are associated with star forma- 
tion episodes triggered by recent merger events. By adopting the 
same procedure used previously to quantify the differences be- 
tween the distributions of the C parameter, we test the differ- 
ences between the D„(4000) distributions shown in Fig. [8] Here 
the results show high confidence in the significance of the dif- 
ferences since none of the 20 randomized realizations exhibit a 
signal as strong as in the data. We conclude that at low metallici- 



Fig. 7. Distribution of the environmental local galaxy density pa- 
rameter and image concentration parameter C for disturbed 
and undisturbed objects in the three metallicity ranges: low 
(solid), medium (dotted), and high (dashed line). 



High Metallicity 



Medium Metallicity 




Fig. 8. Fraction of stellar age indicator (D„(4000)) 
for disturbed normalized to undisturbed galaxies 
(D„(4000)(disturbed)/D„(4000)(undisturbed)) in the three 
metallicity ranges. 



ties, disturbed galaxies have a younger stellar population at least 
at the 95% confidence level. 

We also considered the fraction of disturbed galaxies above 
the median D„(4000) for a given value of M*. By doing so, we 
aim to relate merging events to the relative fraction of young 
stellar populations. In Fig. |9] we show the observed D„(4000) 
vs M, for disturbed and undisturbed galaxies. As can be ap- 
preciated in this figure (see quoted numbers), there is a larger 
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• DISTURBED 
A UNDISTURBED 



% Disturbcd_ 
17.8±2.5% 



10.5 

log(M./M„) 



High Metallicity 



Medium Metallicity 



Low Metallicity 



2 2.5 
(u-r) 



Fig. 9. Distribution of D„(4000) and log(M,/Mo) values for 
galaxies with stellar masses in the range log(M»/M0) ^ 10.0 to 
10.8. The circles correspond to disturbed galaxies and triangles 
to undisturbed galaxies. 



fraction of disturbed galaxies with stellar population ages above 
the median defined by D„(4000) = 0.2125 log(M,/M0) - 0.85 
(shown as a solid line). Thus, our classification of disturbed type 
selects preferentially galaxies with the highest fraction of young 
stellar populations at a given stellar mass content, as expected 
in a merger-induced star formation scenario. We addressed the 
uncertainties in our estimated fractions of disturbed and undis- 
turbed galaxies by moving the median line threshold of 0.3 dex. 
We find that the relative fractions do not change significantly in- 
dicating that the results are not sensitive to the particular choice 
adopted. 

Observati onal analyses of the u - r colour distributions (e.g., 
I Alonso et al.| [2006; De Prop ris et al.l l2005) have reported an ex- 
cess of blue colours in close pair galaxies with signs of inter- 
actions, relative to those measured for galaxies without close 
companions. This finding is associated with both a larger frac- 
tion of actively star-forming galaxies and younger stellar popu- 
lations in disturbed close pair systems. In Fig. [TO] we compare 
the u-r colour distributions of disturbed galaxies with those of 
their undisturbed counterparts. The analysis was performed for 
the three different metallicity ranges previously considered (see 
Table [TJ. We find that in both the high and medium metallicity 
ranges (upper and medium panels, respectively), the colour dis- 
tributions of disturbed objects are more similar to those of undis- 
turbed galaxies than in the lower metallicity range. Interestingly, 
at lower metallicity (lower panel), the colour distributions of dis- 
turbed and undisturbed galaxies exhibit significant differences, 
morphologically disturbed objects having a larger fraction of 
blue colours. We also observe a significant difference in the 
blue tail, where there is a clear excess of blue disturbed galax- 
ies with respect to undisturbed objects. The blue excess can be 
interpreted in terms of substructures associated with externally- 
triggered distortions, indicative of recent interactions or mergers, 
and both enhanced star formation and young stellar populations 
as already shown in Fig. [8] 



Fig. 10. Colour index (u-r) for disturbed (solid) and undisturbed 
galaxies (dashed) in the three adopted metallicity ranges. 



4. Discussion and Conclusions 



Obser ved metallicity gradients in galaxies fe.g.. [ Zaritskv et al 
1994'; Fergus on et al]|1998l; iDutil & Rovlll999t IConsidere et al 



2000) show evidence of less enriched gas in the external re- 
gions. Thus, a large low-metallicity gas reservoir is expected in 
most massive star-forming galaxies. On the other hand, numer- 
ical simulations indicate that galaxy interactions and mergers 
may generate a significant flow of external gas onto the central 
region and trigger star formation. Thus, interactions may be effi- 
cient in lowering the central metal content of galaxies by means 

of an inflow of this external, less enriched gas. 

Similar e ffects have been re ported by Mic hel-Dansac et aTl 

(l2008h and lEUison et al.l (l2008h . where interacting pairs with 
strong morphological disturbance features have lower metallici- 
ties. This is also reinforced by the high ma ss outliers in the MZ 
relation exhibiting distorted morphologies (iPeeples et al ] l2009h . 

The above considerations have led us to analyse the presence 
of possible relics of such process in non-interacting, high-stellar- 
mass galaxies. To this aim, we have performed a morphological 
analysis of high quality images from the MGC of high-stellar- 
mass SDSS galaxies. Objects in this mass range have a more 
clearly defined morphology and so are suitable for exploring 
morphological disturbances associated with merger events. On 
the other hand, the oxygen abundance is a clear indicator of the 
gas phase metallicity and so is sensitive to the presence of re- 
cently formed stars. 

By construction, our subsamples of high, medium, and 
low metallicity objects have similar stellar-mass distributions. 
Moreover, we have confirmed that the subsamples are free of 
strong biases widely different local galaxy density environments, 
or widely different concentration parameters so that we expect 
them to have similar bulge-to-disc ratios. Thus, the subsamples 
analysed comprise galaxies of similar morphology and environ- 
ment. 

Our main result is that we detect a steadily increasing frac- 
tion of morphologically disturbed galaxies, characteristic of be- 
ing merger remnants, with decreasing O/H. The three bins anal- 
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ysed consist approximatively of 15%, 20%, and 50% disturbed 
galaxies (high, medium, and low metallicity, respectively) 

Moreover, the ratio of disturbed to undisturbed distributions 
of the D„(4000) parameter in both the medium and high metallic- 
ity bins are similar, the D„(4000) values corresponding to a rel- 
atively old stellar population in disturbed objects. This contrasts 
with the results for the low metallicity bin where disturbed ob- 
jects have a substantially younger stellar populations than their 
undisturbed counterparts. In addition, the colours distributions 
exhibit similar trends, showing that low metallicity galaxies with 
a disturbed morphology are bluer than those that are undisturbed. 

The bluer colours and lower D„(4000) values in low metal- 
licity, morphologically disturbed objects, suggest externally- 
triggered distortions driven by a recent interaction or merger, 
with an associated enhanced star-formation rate and a predomi- 
na nce of a you ng stellar population. 

iRupke et al.l (12008) showed that ultra-luminous infrared 
galaxies (ULIRGs) have systematically lower O/H values than 
star-forming galaxies of similar stellar mass content, consistent 
with a merger-driven inflow of low metallicity gas. Our results 
show, in addition, that this mechanism is not only present in 
ULIRGs but also in high stellar mass galaxies with morphologi- 
cal evidence of a merger event. 

The most natural explanation of our results is that merger- 
induced low metallicity gas inflows have occurred from the ex- 
ternal regions of high stellar mass galaxies. Our results also sug- 
gest that the timescale of the chemical evolution (the combined 
effects of the enrichment of the interstellar medium by SNII and 
these inflows) is at least as long as the timescale of the tidally 
induced morphological disturbances. Thus, galaxies with a dis- 
turbed appearence may contain young stellar populations and yet 
low gas-phase abundances in their central parts. Our results sug- 
gest that, at least in part, the scatter in the MZR at high stellar 
mass is caused by interactions and mergers. 
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